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ABSTRACT 

We study the galaxy number density (smoothed on a 5/i -1 Mpc scale) around 
18 low-redshift Lyman-alpha absorbers previously observed with HST. The ab- 
sorbers lie in the foregrounds of Mrk 335, Mrk 421, Mrk 501, I Zw 1, and 3C 273, 
all within regions where there are now complete redshift surveys to mz w — 15.5. 
We construct a smoothed galaxy number density field from the redshift survey 
data and determine the distribution of densities at the Lyman-alpha absorber 
locations. We also find the distribution of galaxy number density for a variety 
of test samples: all galaxy locations within the Center for Astrophysics Redshift 
Survey (CfA2), CfA2 galaxy locations along randomly selected lines of sight, and 
randomly chosen redshifts along random lines of sight. 

The Lyman-alpha absorbers are present in dense regions of the survey, but 
occur far more frequently in underdense regions than do typical luminous galax- 
ies. The distribution of smoothed galaxy density around the Lyman-alpha ab- 
sorbers is inconsistent at the 4<r level with the density distribution around survey 
galaxies. It is highly consistent with a density distribution at randomly chosen 
redshifts along random lines of sight. This supports earlier evidence that the 
nearby, low column density (logA^ji ^ 14) Lyman-alpha forest systems are spa- 
tially distributed at random; they are not well correlated with the local large-scale 
structure. 

Subject headings: large-scale structure of the universe - - quasars : 
absorption lines — intergalactic medium 



1. Introduction 

The spectra of high-redshift QSOs display a dense network of absorption 
lines blueward of the QSO Lyct emission, first recognized by Lynds (1971) as low 
column density Lya absorption in the rest frame of intervening discrete systems 
(cf. reviews by [Bajtlik 1993] ; |Weymann 19931) . With the advent of the HST, this 
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"Lya forest" could be studied at wavelengths below the atmospheric UV cutoff, 
corresponding to absorbers with redshift z <^ 1.6. Because the Lya forest reveals 
a considerable population of very tenuous structures (N Hl > 10 12 cm -2 ) across 
cosmological distances (z ^ 4), there has been great interest in determining the 
nature of the absorbers: galactic, intergalactic, or both. 

Early studies of the redshift distribution of Lya forest clouds at high z re- 
vealed that the metal-poor population was much less clustered than the present- 
day galaxies QSargent et al. 1980|) . Furthermore, there were suggestions that 
the Lya forest might be an entirely separate population from the absorbers at 
larger column density (Nm > 10 17 cm -2 ), which have associated metal-line ab- 
sorption and cluster like galaxies ( [Sargent, Boksenberg, fc Steidel 1988| [lytler 



1987[). Proposed explanations of a non-galactic Lya forest include pressure- 



confined clouds in the IGM ( |Sargent et al. 1980| ; |lkeuchi fc Turner 1991| ) and 
CDM minihalos fllkeuchi 19851 |Rees 19861 ). Recent results from the HIRES 
spectrograph on the Keck 10m telescope still admit the possibility of a dual- 
population (galactic and intergalactic) Lya forest. Highly resolved metal-line 
absorption in high-z damped-Lya systems appears to originate in thick, rapidly 
rotating (^ 200 km s -1 ) protogalactic disks ( Prochaska fc Wolfe 1997| ). The 
mean metallicity of high-z Lya clouds declines sharply at low column densities, 
from [C/H]~ -2.5 at N m > 10 14 - 5 cm" 2 to an upper limit of [C/H]< -3.5 for 
10 13 - 5 < N m < 10 14 cm" 2 ( [Lu et al. l99g| ). 

With the HST discovery of several z < 0.16 Lya clouds along the 3C 273 
sightline ( [Morris et al. 199"T| ; [Bahcall et al. 199~T| ), it was now possible with pencil- 
beam redshift surveys to search for far less luminous galaxies coincident in redshift 
space with the low column density absorbers. Morris et al. (1993) completed a 
B ~ 19 survey of ~ 4 square degrees around 3C 273 that detected no galaxies 
within 1.5/t Mpc of any absorber. They concluded that the absorption is 
extragalactic in origin. A further deep imaging survey of 53' x 52' around 3C 273 
by Rauch, Weymann, & Morris (1996) to a limiting central surface brightness of 
fi r = 26.4 revealed no low surface brightness galaxies within projected separations 
of 80/i" 1 kpc and 123/iT 1 kpc, respectively, at the two nearest absorber redshifts. 

Lanzetta et al. (1995), on the basis of an r < 21.5 survey within 1'3 of several 
HST/FOS target QSOs ([Bahcall et al. 1993|) , concluded instead that the Lya 
forest is largely due to extended galactic gaseous halos QBahcall fc Spitzer 1969| ) 
or gaseous disks ( Maloney 1993 ; Hoffman et al. 1993j ) of size w lQOh^ 1 kpc and 
roughly unit covering factor. Their survey spanned a much greater cumulative 
pathlength, Az m 2.86, than the Az ~ 0.16 sampled by the Morris et al. (1993) 
study of the 3C 273 foreground. Lanzetta et al. detected 30 galaxies at impact 
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parameters <^ 160/i _1 kpc from the QSO lines of sight (LOS), closer than any 
detected galaxy to the 3C 273 LOS. They found: (1) a strong increase in the 
probability of a galaxy-absorber association for galaxy impact parameters under 
\60h~ 1 kpc. They also noted a statistical ant i- correlation between the equivalent 
width of Lya absorption and the impact parameter of a galaxy-absorber pair. 
They argue that at z <^ 1, the fraction of Lya absorption systems arising in 
galactic halos is at least 0.35 ± 0.10, and probably 0.65 ± 0.18. It is notewor- 
thy that their HST/FOS Lya absorber sample was limited to equivalent widths 
W^ya ^ 0.3A, larger than almost all the ffST/GHRS detections and probably not 
representative of the extremely metal-poor population, with Nhi < 10 14 cm~ 2 . 
Mo and Morris (1994) used Monte Carlo simulations to show that the galactic- 
halo fraction of 3C 273 absorbers, all with ^ 10 14 cm -2 , could only be ~ 20% 
given the observed number of absorbers and the weak correlation between ab- 
sorbers and galaxies on large scales. 

In a different approach, Shull, Stocke, & Penton (1996, hereafter SSP96; 
Stocke, Shull, fc Penton 1995|) obtained HST/GHRS spectra of four nearby 



bright AGNs within the boundaries of existing wide-angle redshift surveys to 
determine whether Lya absorbers nearby are correlated with large-scale struc- 
ture, as first suggested by Oort (1981). They detected a total of 11 Lya ab- 
sorbers along these sightlines within a GHRS wavelength range corresponding to 
1500 ^ cz <^ 10500 km s . Using the wavelet analysis of Slezak, de Lapparent, 
& Bijaoui (1993) to subdivide their AGN sightlines into "void" and "superclus- 
ter" regions, SSP96 found four of their eleven absorbers in voids. Because their 
sightlines probed an approximately equal pathlength of void and supercluster, 
binomial statistics gave a fair probability (19%) that the z < 0.035 absorbers 
are uniformly distributed and not biased towards superclusters. From the rela- 
tion between Lya equivalent widths and nearest-galaxy distance, SSP96 argue 
that the W\jy a < O.lA absorbers are distributed in a manner statistically indis- 
tinguishable from clouds randomly placed with respect to galaxies. Although 
the nearest-galaxy distances of their Wi, ya > O.lA subsample are broadly consis- 
tent with the Lanzetta et al. (1995) trend extended to larger impact parameters 
(^ 500/i -1 kpc), it is unlikely that the Lya absorbers could be physically associ- 
ated with galaxies at that separation. 

We use the large-scale smoothed galaxy number density to compare the dis- 
tribution of z < 0.035 Lyct absorbers with that of typical luminous galaxies. 
Whereas SSP96 identified the absorbers as either "void" or "supercluster", we 
quantify the large-scale structure in the vicinity of the absorbers. The smoothed 
galaxy density field is a more robust indicator of absorber environment than the 
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nearest-galaxy distance, particularly on scales larger than 1 Mpc. 

The technique we develop requires a wide-angle redshift survey overlapping 
the AGN sightlines of interest and adequately sampling the galaxy population at 
the Lya absorber redshifts. Thus we restrict ourselves to the cz < 10500 km s _1 
absorbers along the sightlines to 3C 273 and the four AGNs of SSP96. For the 
18 Lya absorbers meeting these specifications, we use the galaxy number density 
field from the CfA2 Redshift Survey flGeller fc Huchra 19891 ; IVogeley et al7T994l 
Marzke, Huchra, fc Geller 1994| , hereafter MHG94) and adjacent regions also 
surveyed to m Zw = 15.5 flMarzke et al. 1996| ; |Grogin, Geller, fc Huchra 1998|) . 
Even with this modest sample of absorbers, we obtain rather tight constraints 
on the relation between the local Lya absorbers and the large-scale structure 
marked by intrinsically luminous galaxies. 

Section describes the nearby Lya absorption systems and the galaxy redshift 
surveys surrounding them. Section |3] describes the density estimation technique 
and shows maps of the galaxy density field toward the five target AGNs. We also 
derive the cumulative distribution function (CDF) of the densities around the 18 
Lya absorption systems. In §(| we investigate the implications of the Lya density 
CDF by comparing it with CDFs for several test populations. We also compute 
likelihood intervals for models of the absorbers' spatial distribution. We address 
sources of uncertainty in §|5|, and estimate the improved constraints likely with a 
larger sample of nearby Lya absorbers. We conclude in §0. 



2. Data 

2.1. Ly man-alpha Absorption Systems 

Our sample of local Lya absorbers comprises seven systems with cz < 10500 



km s 1 detected towards 3C 273 ( [Morris et al. 1993]) and eleven systems detected 



in the sightlines toward the AGN Mrk 335, Mrk 421, Mrk 501, and I Zw 1 
( (5SP96| ). The 3C 273 observations extend down to 1215A ([Weymann et alT 



1995|) and include Lya forest systems well into the damping wing of Galactic 



Lya absorption. The AGN observations span a more limited wavelength range of 
1222A-1259A, corresponding to a redshift window of 1500 ^ cz <^ 10500 km s _1 
for the detection of Lya absorption systems. Three of the four AGN have cz < 
10500 km s _1 , further reducing the path length probed by these sightlines. Table 
[l] summarizes the relevant data for each of the Lya sources, including the source 
redshift and the redshift range for absorber detection. In Table |2] we list the 
redshifts and equivalent widths for the 18 Lya absorbers. 
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2.2. Redshift Surveys 

The CfA2 Redshift Survey (|Geller fc Huchra 1989|) provides the broad sky 
coverage and dense sampling necessary to evaluate the galaxy number density, 
smoothed on scales <^ 5/i _1 Mpc, around Lya absorbers with cz <^ 10000 km s _1 . 
Although there have been much deeper, pencil-beam surveys along QSO sight- 
lines ( fYlorris et al. 1993] ; [Lanzettaet al. 1995| ) to link Lya absorption with specific 
galaxies, these subtend much too small an angle to evaluate nearby densities on 
5/i _1 Mpc scales. SSP96 chose their four AGN partly because they are within the 
CfA2 boundaries as analyzed by Vogeley et al. (1994) and MHG94: 8 h < a < 17 h , 
8?5 < 5 < 44?5 (CfA2 North) and 20 h < a < 4 h , -2°. 5 < 5 < 42° (CfA2 South). 
Redshifts around these four AGN sightlines are contained in Huchra et al. (1990), 
Huchra, Geller, & Corwin (1995), Huchra, Vogeley, & Geller (1998), and are 
available through NEDQ. 

We have also surveyed a 2 h x 12° region nearly centered on 3C 273 ( |Grogin 



Geller, fc Huchra 1998; ). With « 1060 redshifts in the region ll h 30 m < a < 
13 h 30 m and -3?5 < 5 < 8?5, this survey extension is 99.5% complete to the CfA2 
magnitude limit (mz w = 15.5) and 98.1% complete to the Zwicky catalog limit 
(mzw = 15.7). This additional survey yields estimates of the global density at the 
locations of the 3C 273 Lya absorbers. In addition to the 3C 273 field redshifts, 
we also incorporate ~ 1700 Zwicky galaxy redshifts at low Galactic latitude from 
the Galactic Plane Survey by Marzke, Geller, & Huchra (1996). These additional 
redshifts outside the CfA2 boundaries prevent substantial underestimation of the 
galaxy density at locations within a smoothing length of the survey boundaries. 



3. Estimating the Smoothed Galaxy Number Density 

Our goal is to compare the distribution of Lya absorbers with the large- 
scale distribution of galaxies. If the absorbers sample the large-scale structure 
marked by luminous galaxies in a redshift survey, the distribution of smoothed 
galaxy number densities around the absorbers should be consistent with the 
distribution of densities around survey galaxies. Conversely, if the Lya absorbers 
arise in a host population distributed independently of the large-scale structure, 
the densities around the absorbers should resemble the densities around similarly 
distributed locations in the survey volume. 



1 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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In § 1371] we describe the smoothed galaxy number density field constructed 
from the redshift survey data of §2.2j . We show maps of the density field towards 
five AGN with low-redshift Lya absorbers in §[J]2]. Finally we combine the den- 
sities at each of the 18 absorber locations into a cumulative distribution function 



§3.3) for comparison with density CDFs of various test populations in 



3.1. Method 

To evaluate the galaxy density around the local Lya absorbers, we first trans- 
form the point distribution of the CfA2 Survey into a continously defined number 
density field throughout redshift space. We smooth each galaxy in redshift space 
by a unit-normalized Gaussian kernel W of width a = 5/iT 1 Mpc: 

W(x- x gal ) = (2na)- 3/2 exp (|x - x gal | 2 /2a 2 ) . (3-1) 

We choose a 5/i -1 Mpc smoothing length to coincide with the galaxy-galaxy 
correlation length ([Peebles 1993| ; Marzke et al. 19951 ; |Jing, Mo, fc Boerner 1998[ ) 



and with the pairwise velocity dispersion in the survey flMarzke et al. 1995|) . In 



§ p72| we discuss the sensitivity of our results to this choice. 

We make no attempt to remove peculiar velocity distortions (cluster "fingers" , 
etc.) from the redshift survey or from the Lya absorber velocities. In principle, 
they should share the same local velocity field. We do introduce the standard 
heliocentric to Galactocentric correction, 

cz = czq + (300 km s -1 ) sin I cos b, (3-2) 

for an object at Galactic longitude I and latitude b. We place each object at a 
comoving distance r appropriate for a go — 0.5 universe with pure Hubble flow: 

r(z) = (^)[l-(l + z)-^]. (3-3) 

We thus underestimate spatial overdensities associated with clusters, which are 
broadened in the radial direction. Our smoothing kernel effectively washes out 
peculiar velocities ^ 500 km s _1 , close to the 540 ± 180 km s _1 pairwise velocity 
dispersion measured by Marzke et al. (1995) for the combined CfA2 and SSRS2 
( |da Costa et al. 1994|) surveys. 

Because the CfA2 Redshift Survey is flux-limited, an increasing fraction of the 
galaxies at larger redshift fall below the magnitude limit and do not appear in the 
survey. In computing the density field, we compensate for the magnitude-limited 
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sample by assigning each galaxy a weight l/if>, where the selection function ip is 



■0(o!, 5, z) 



4>{M) dM 



0(M) dM 



(3-4) 



Here M\ [m is the effective absolute magnitude limit at the galaxy position and 
4>{M) is the differential luminosity function (LF). For Mn m > M cut , we assign 
galaxies unit weight. All magnitudes are Zwicky magnitudes. Numerical values 
for absolute magnitudes implicitly include the /i-dependence in equation ( |3-3[ ) . 

MHG94 fit the CfA2 LF to a Schechter function SF ( fcchechter 1976| ), con- 
volved with a Gaussian error of om = 0.35 mag ( [Huchra 1976] ) in the Zwicky 
magnitudes: 

SF (M) = ^(0.41nl0)10 a4(M *- M)(1+Q) exp[-10°- 4(M *- M) 
1 



0CfA2(M) 



2na M J - 



S F(M')exp — (M' — M) /2a M ) dM'. (3-5) 



We adopt the values 0* = 0.04 (Mpc//i)~ 3 , M* = -18.8, and a = -1.0 found 
by MHG94 using M cut = —16.5. We discuss the sensitivity of our results to 
the LF parameters in § |5.1| . For computational convenience in determining ip, we 
replace the convolution of equation (|3-5|) with (ficfA2(M) w 0sf(-^ + 0.1 mag). 
This approximation recovers the true ip to better than 5% for cz ^ 12000 km s _1 . 

For a galaxy at position (a, S) and at luminosity distance D L (z) = (1+z) r(z), 
we estimate Mi im according to 



M lim (a, 6, z) = m hm - 5 log 



1 + z) r(z) 
l/t^Mpc 



25 - Am K (z) + Am ext (a, 5). (3-6) 



In equation (|3-6| ), mii m is the CfA2 flux limit, Amx is a i^-correction, and Am ex t 
is a correction for Galactic extinction. Photoelectric photometry of Zwicky galax- 
ies ( [lakamiya, Kron, fc Kron 1995| ) suggests that Volume I of the Zwicky catalog 
( [Zwicky et al. 1961-1968p goes ^0.4 mag fainter than the other volumes at the 
CfA2 magnitude limit, mz w = 15.5. To correct for this Volume I scale error, we 
adopt m\ im = 15.9 for all galaxies with 5 < 14?5 in CfA2 North and the 3C 273 
region extension. 

Lacking morphological types for the majority of CfA2, we apply a generic K- 
correction, Am K (z) = 3z, appropriate for type Sab (|Pence 1976|) . Our error in 
the K-correction will be small due to the predominantly low redshifts (z <^ 0.05) 
in the survey. To obtain our correction Am ex t(tt, 5) for Galactic extinction along 
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a particular line of sight, we first interpolate the Hi map of Stark et al. (1992). 
We then convert from Hi to reddening with the relation (N(H.i)/E(B — V)) = 
4.8 x 10 21 cm" 2 mag -1 ( ^ombeck 1990[ ), and adopt a standard extinction law 
Am cxt = A B = A.0E{B-V). 

We compute the smoothed galaxy number density n at a given point r = 
(a, 5, r(z)) by summing the contributions from all i galaxies in the survey: 

W(r - n) 



n(r) = 



E 



fa 

W(r-n) / 4>(M)dM 



<f>(M) dM 



<j>{M) dM 

For the CfA2 Survey, MHG94 derive a mean density n = J^ ut 0(M) dM = 
0.07 (Mpc//i)~ 3 with M cu t — —16.5. This analysis yields a distribution of dimen- 
sionless galaxy density contrasts, n(r;)/n, at the absorber positions r$ which may 
be compared with the n/n distribution of various control samples. 



3.2. CfA2 Density Maps toward Local Lyman-alpha Systems 

For each of the five sources in Table [l], we display the results of our density 
mapping in two forms . The first representation is similar to the standard "wedge 
plot" format for redshift survey data: we project galaxy positions over a narrow 
declination range onto a circular section with R.A. as the azimuth coordinate and 
redshift as the radial coordinate. We add the density contrast as a superposed 
isodensity contour map, made by sampling the density at regular 3/i -1 Mpc 
intervals on the midplane (or, more accurately, the mid-cone) of the wedge plot. 

We set the outer boundary of the wedge plots at 15000 km s" 1 , well into the 
region where shot noise dominates the density estimate. We choose a "thick- 
ness" of six degrees in declination, centered on each LOS. Our 5/i _1 Mpc den- 
sity smoothing length subtends this angle at a redshift of ~ 5000 km s _1 . For 
cz ^ 5000 km s _1 , there is density variation with declination that is lost in pro- 
jection. To showcase the large-scale structure toward each target, the wedge plots 
have a 30° opening angle centered on each LOS, corresponding to (2 h / cos<5) of 
right ascension at the central declination S. Our density smoothing length sub- 
tends this angle at a redshift of (1000/ cos 5) km s _1 . We accommodate the 
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dynamic range of the density map with underdensity contours (dotted) in linear 
decrements of 0.2n and overdensity contours (solid) in logarithmic increments 
corresponding to n, 2n, An, 8n, etc. In addition to the surveyed galaxies (small 
squares), we also plot the locations of the detected Lya absorbers (large trian- 
gles) along the line of sight (dashed) to the background source (large circle, if 
cz src < 15000 km s" 1 ). 

For the second display, we sample the density at 100 km s _1 radial intervals for 
cz < 15000 km s -1 along the line of sight to the selected background source. We 
then plot the density profile (n/n) along the LOS (dashed curve) versus redshift. 
We indicate the redshifts of the detected Lya absorbers (dotted lines) as well 
as the redshift interval over which the absorbers could have been detected (solid 
lines). If the upper redshift limit coincides with the source of Lya emission, we 
indicate this coincidence with a doubled solid line. Within this redshift interval, 
we also show the ±lcx uncertainty in the local density estimate (solid curves). 
This density error estimate is the quadrature sum of the fractional errors due to 
sampling variance (cf. § |5.3| ) and the uncertain survey LF (cf. § |5.1|) . 



3.2.1. 3C 273 



For the wedge plot of the 3C 273 region (Figure [I]), we center the right as- 
cension on the 11^5-1 3^5 range of the CfA2 survey extension (|Grogin, Geller, 



fc Huchra 19981) instead of the ±15° about 3C 273 (a = 12 h 26 m 6). Because 
the survey extension is complete to the Zwicky catalog limit, mz w = 15.7, we 
also show the locations of the fainter galaxies (small crosses) to compare with 
the Lya systems. We use only the brighter mz w < 15.5 galaxies to derive the 
density estimate. The uncertainty in the density estimator rises dramatically 
for cz ^ 12000 km s _1 (cf. |5]3|), and exceeds ±1 dex at the two absorption sys- 
tems near 15000 km s" 1 . We therefore only consider the 3C 273 absorbers with 
cz < 10500 km s _1 , coincident with the velocity maximum for the other AGN 
lines of sight. 

The seven 3C 273 absorbers with cz < 10500 km s -1 inhabit a wide range of 
environments. The two nearest Lyct systems are near the Virgo cluster, and thus 
are difficult to identify in the crowded wedge plot. Their redshifts are clearly 
indicated in the associated density profile. These two systems inhabit some of 
the highest-density regions in the LOS. In contrast, two of the other absorbers 
inhabit the low-density extreme of the LOS, a large void at cz pa 9500 km s _1 . 
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3.2.2. Mrk 335 

Figure shows the Mrk 335 region. The local galaxy densities at the Lya 
absorber locations span the density range in the LOS redshift window, including 
a pair with cz m 2350 km s _1 near the density minimum. 

3.2.3. Mrk 421 

Figure |3] shows the Mrk 421 region. This LOS is intriguing in that only one 
absorber appears within the ~ 7500 km s" 1 sensitivity window; it resides in a 
relatively under dense environment. 

3.2.4. Mrk 501 

The proximity of Mrk 501 (a = 16 h 52 m 2) to the CfA2 North boundary at 17 h 
poses an added challenge to determination of the densities surrounding its Lya 
absorbers. Fortunately there is an essentially complete redshift survey ( [Marzke et| 
al. 199(j| ) beyond this boundary, mitigating edge effects on the density estimate. 
In the region 17 h < a < 18 h and +34° < 5 < +46°, there are redshifts for 140 of 
the 150 Zwicky galaxies with mz w < 15.5. 

Rather than cropping the Mrk 501 wedge plot at the CfA2 survey boundary, 
we set the upper R.A. limit in Figure ^ at 18 h in order to show the large-scale 
structure east of the LOS. As noted above, there are redshifts for more than 
90% of the galaxies in the plotted region with a > 17 h . We therefore have high 
confidence in the density estimates along the Mrk 501 LOS. 

Interestingly, there are no absorbers in the region of large overdensity beyond 
cz ~ 8000 km s" 1 associated with the Hercules supercluster. On the contrary, 
two of the absorbers are near the LOS density minimum, in a foreground void to 
the supercluster. A third is on the outskirts of this void, in a region of middling 
density. 

3.2.5. I Zw 1 

Figure shows the I Zw 1 region. The two nearest absorption systems are 
located in a relatively shallow void; the third system appears to be associated 
with the galaxy wall behind it at cz « 5500 km s" 1 . Beyond this wall, the LOS 
passes through a larger and emptier void containing no detected Lya systems. 
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3.3. Density Distribution for the Lyman-Alpha Absorber Sample 

Figure || shows the cumulative distribution function (CDF) of the densities 
surrounding the 18 Lya absorbers (solid). We plot the density contrast (n/n) 
on a logarithmic axis. Error in the overall normalization of (n/n) translates 
all CDFs along the x-axis but does not change their shape. The distribution of 
(n/n) for the Lya absorbers is broad, with fully half the absorbers in substantially 
underdense regions of the survey: (n/n) < 0.4. This considerable fraction of low- 
(n/n) absorption systems is the principal discriminant among the density CDFs 
of test samples in §f|. 

4. Interpretation 

To understand the physical implications of the galaxy density distribution 
around the sample of nearby Lyct absorption systems, we define several test pop- 
ulations. We then calculate the local densities surrounding each of the locations 
in the test samples. Finally, we employ a two-population Kolmogorov-Smirnov 
(K-S) statistic to compare the galaxy number density distribution of the Lya ab- 
sorbers with the density distributions from the test samples of: (1) CfA2 galaxies 
(§|4.1|, § [4.2|) and (2) randomly selected locations within the redshift survey bound- 
aries (§ |4.3|) . To construct these test samples, we must also simulate the radial 
sampling bias caused by the redshift windows of the Lya forest observations 
(Table 

4.1. Comparison with All CfA2 Galaxies 

For this test sample, we determine the density contrast (n/n) at the locations 
of each CfA2 galaxy (within the boundaries listed in §[T2] and including the 3C 
273-region extension). To account for the radial sampling bias of the Lya sight- 
lines, we assign each of the densities in the sample a weight W equal to the number 
of redshift intervals in Table p] which contain the redshift of the corresponding 
galaxy. The value of W therefore varies between zero (cz > 10660 km s _1 ) and 
five (1713 < cz < 7901 km s" 1 ). 

When using CfA2 galaxies to trace the density field, the survey magnitude 
limit causes an additional radial sampling bias. Even though we have corrected 
for the survey selection function ip in determining (n/n), the CfA2 galaxies at 
larger redshifts are sampling (n/n) much more sparsely because we see propor- 
tionately fewer of them per unit volume. We attempt to correct for this sampling 
bias by assigning (n/n) for each galaxy an additional weighting factor of 
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Thus each (n/n) in the sample has weight (W/i/j). 

The weighting is not an entirely satisfactory solution to the radial bias prob- 
lem, because the unseen fraction (1 — ip) of galaxies at higher redshift are not all 
located at the positions of the observed fraction ip. The fixed smoothing length 
of the density estimator results in pronounced local maxima in the density field 
at the locations of higher-redshift galaxies in the survey. We therefore tend to 
overestimate the high-density end of the (n/n) distribution when assuming that 
a particular galaxy's (n/n) is representative of its (ip -1 — 1) unseen neighbors. 

Figure |^ shows the cumulative distribution function (CDF) of (n/n) for all 
CfA2 galaxies, with (W/ip)- weighting (dotted) as well as with W- weighting alone 
(dashed). For comparison, we also plot the density CDF for the Lya absorbers 
(solid). Clearly the Lya absorbers occur far more frequently in underdense re- 
gions. The K-S probability that the Lya CDF and the (jy/?/>)-weighted galaxy 
CDF are drawn from the same underlying distribution is only 3 x 10~ 8 , a 5.1a 
rejection of the null hypothesis. The K-S probability for the W-weighted CDF is 
9 x 1CT 4 , a 3.3a rejection. 

4.2. Comparison with CfA2 Galaxies in Random LOS Cylinders 

The HST observations of local Lya absorbers sample the population along 
narrow sightlines through redshift space. This geometric sampling is not repro- 
duced in the volume-filling, complete CfA2 galaxy sample ( §4.1| ). Discrete lines 
of sight through the survey are less likely to sample rich galaxy clusters because 
the clusters within ~ lOO/i -1 Mpc cover only a small fraction of the sky. We 
therefore define another test population by selecting random LOS through the 
CfA2 survey and extracting the galaxies within cylinders of diameter 5h~ l Mpc 
(the density smoothing length) around those LOS. 

For each of the five actual AGN sightlines (Table |l|), we assign the corre- 
sponding redshift interval to one fifth of the randomly generated LOS cylinders. 
To obtain a suitably large test sample of 5000 galaxies, we use ~ 150 LOS. We 
again correct for the selection function sampling bias (cf. § f4.1| ) by assigning each 
galaxy's (n/n) value a i\)~ x weight in the census. 

Figure [7] shows the CDF of (n/n) for the random-LOS CfA2 galaxies, 
weighted (dotted) as well as unit-weighted (short dash). We again see a disparity 
between the Lyct and CfA2 galaxy CDFs, though perhaps less pronounced than 
for the full-CfA2 sample. The sharp rise in the unit-weighted CDF at (n/n) m 4 
is caused by a combination of: (1) a preponderance of Virgo cluster galaxies, 
all at (n/n) « 4, appearing in the one-fifth of LOS cylinders which extend to 
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z = and (2) a lack of (n/n) > 4 galaxies because of the low probability that a 
5/i _1 Mpc cylinder intersects a CfA2 cluster (other than Virgo), or even a galaxy 
with cz ^ 8000 km s _1 . This feature is much less apparent after we apply the 
selection-function weighting correction. 

The K-S probability that the Lya CDF and the ?/>~ ^weighted galaxy CDF 
are drawn from the same underlying distribution is 3 x 10 -5 , a 4.2<r rejection. 
The K-S probability for the unit-weighted CDF is 2 x 10~ 4 , a 3.7a rejection. 



4.3. Comparison with Randomly Generated Velocities along 

Random LOS 

We next construct test populations of Lyct absorption systems with randomly 
chosen redshift-space coordinates along randomly selected LOS. We first adopt 
a contstant linear density N Lya (cz) = 1/(2500 km s -1 ) for our mock absorbers, 
roughly equivalent to the observed mean incidence for the five selected HST lines 
of sight (18 absorbers over ~ 42000 km s _1 ). Because all these sightlines have 
2 max < 0.04 and span a redshift range Az « 0.03, we ignore the linear density 
gradient dN-^ ya (z)/dz. As in the previous case, we assign radial intervals to these 
random LOS by selecting evenly from the five redshift windows of Table [|. 



4-3.1. Single-population Model 

For a given random LOS through CfA2, we sample the smoothed galaxy den- 
sity at lh^ 1 Mpc radial intervals to obtain a density- vs.-redshift profile n(z). We 
take this density profile, raised to an exponent f3, as the probability distribution 
function for the mock-Lya redshifts: 

P{z) dz = ,, max ' • (4-1) 

nP{z') dz' 



By drawing a sample of mock absorbers from a (3 = (or P(z) = const.) model, 
we can test the hypothesis that low column-density Lya absorbers are distributed 
uniformly throughout space and completely uncorrelated with local galaxy den- 
sity. Similarly, we could test the hypothesis that the Lya systems reside in density 
environments typical of CfA2 galaxies by drawing a sample of mock absorbers 
from a /3 — 1 model. 

The distribution power-law index /3 may be viewed very roughly as the com- 
bination of a cloud-creation index (3 C and a cloud-destruction index /3 'a- When 
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mergers, etc., destroy low column- density absorbers in regions of high galaxy 
density (/3d > 0), we might expect an anticorrelation of absorber location with 
local galaxy density ((3 < 0). To estimate (3 in a merger-destruction scenario, we 
note that the rate of cloud-galaxy interactions scales as the product of the local 
galaxy cross section (oc n) and the peculiar velocity of the clouds through the 
redshift space (oc n^ 1 / 2 ). If the clouds were produced in direct proportion to the 
local galaxy density ((3 C = 1), then j3 d w 3/2 would imply (3 ~ —1/2. 

Only in the case (3 = do we explicitly simulate a sample of ~ 5000 mock 
absorbers to compare with the observed set of Lya absorbers. Figure [j] shows 
the density CDF of the 5000 randomly located absorbers (long dash), which may 
be compared with the CDF for the observed Lya systems (solid). The agreement 
between the two CDFs is remarkable, particularly when compared with the CfA 
galaxy CDFs in Figs. | and [7[ The K-S probability that the observed Lya 
densities share the same underlying distribution as the (3 = model is 0.20, 
consistent at the 1.3<r level. This probability only varies by ±0.01 for multiple 
realizations of the randomly generated sample, indicating that a 5000-member 
sample is sufficiently large for our analysis. 

We check the sensitivity of this result to the adopted linear density of ran- 
domly generated absorbers. When we vary the linear density by ±30%, the K-S 
probability varies by ±0.03, larger than the sampling variance (0.01) but still 
statistically insignificant. We also test a simulation in which the redshift interval 
and the linear density of absorbers along a random sightline are set to the ob- 
served values for one of the five AGN sightlines. This sampling compensates for 
varying sensitivity for Lya detection between sightlines. The analyses of the 3C 
273 HST spectra yielded detections of lower equivalent-width absorption (^ 30 
mA) features than SSP96 detected for their AGN. Not surprisingly the 3C 273 
LOS has a significantly higher linear density of absorbers (cf. Tables [l] and ||]). 
The K-S probability for this variable-linear-density model is the same as the 0.20 
probability for the fixed-density model to within the sampling variance. The 
match between the density CDFs of the observed Lya clouds and the randomly 
chosen locations in redshift space is robust against reasonable variation in the 
adopted linear density of mock absorbers. 

We next use a maximum-likelihood technique to estimate the range of j3 con- 
sistent with the observed density distribution of local Lya systems. If equation 
( |1- 1| ) describes the actual redshift distribution of local Lya absorbers, then the 
likelihood of seeing the observed redshift distribution (Table |2]) is the product of 
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the individual probabilities for each absorber at Vf 



m = n = n ^ / — ( 4 - 2 ) 

u(a i} 5 i} z') dz 1 



Given the prior assumption that the underlying distribution of Lya locations 
follows a /3-model, then the quantity — 21n[£(/3)] should be distributed as x 2 m 
(3. 

In Figure |], the solid curve shows the model likelihood for the range —1.5 < 
(3 < 1.5. We obtain a maximum likelihood value of j3 = —0.02 ± 0.23 from the 
18 nearby Lya systems. A single-population Lya forest that traces the local 
large-scale structure (j3 = 1) is inconsistent with our result at the 4<r level. This 
result may be compared with the ~ 4cr inconsistencies between the density CDFs 



of the Lya absorbers and the galaxy test populations (§§[44], fL2|) , which should 
have a (3 — 1 distribution throughout the survey. 



4-3.2. Dual-population Model 

Mo and Morris (1994) speculated that the Lya absorbers may include two 
distinct populations: one population associated with the gaseous extended ha- 
los of galaxies, and another population of isolated, unclustered clouds. In this 
model, we might expect the galactic-halo fraction / gal to have local densities con- 
sistent with the densities around typical galaxies in the survey {(3 — 1) and the 
unclustered fraction to have local densities consistent with randomly chosen lo- 



cations in the survey ((3 = 0). Equation ( |4-1|) then yields a distribution function 
of dual-population absorbers of the form: 

P(z) dz = ^1 <Z) ^ + / (Wsal)tfe . (4-3) 



n(z') dz' 

In analogy with equation ( |4-2j) we can define a likelihood function for the observed 
clouds to be drawn from a model with galactic fraction / ga i: 

£(/ g al)=n P ( r -/ g al), (4-4) 



where P is now given by equation ( |4-3| ). 

In Figure || the solid curve shows the likelihood of a dual population model 
with fraction / ga i distributed like the CfA2 galaxies ((3 = 1) and the remainder 
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distributed uniformly ((3 = 0). The maximum-likelihood value for the galactic- 
type fraction is / ga j = 0.00, with a la upper limit of / ga i = 0.24. This result is 
inconsistent at a 99% (2.6a) confidence level with the upper estimate of / ga i ~ 
0.65 given by Lanzetta et al. (1996). Their lower limit of / ga i ~ 0.38 is only 
excluded at a confidence level of 88% (1.5a). 



5. Discussion 

Our Lya local density analysis suffers from a number of sources of uncertainty 
which increase the error bars in our results of §^j. Because we assign each galaxy 
a weight derived from the CfA2 selection function, the density estimator is sen- 
sitive to uncertainty in the shape of the CfA2 luminosity function. The density 
estimator, constructed from a point distribution of galaxies, is subject to increas- 
ing shot noise at larger redshift. The shot noise dominates at cz ^ 10000 km s -1 , 
where the mean distance between survey galaxies rapidly exceeds our 5/i _1 Mpc 
smoothing length. For this reason, we may worry that the parameter constraints 
are also sensitive to the adopted smoothing length. Here we quantify some of 
these sources of uncertainty and estimate the improvements in Lya distribution 
constraints which would result from a larger sample of nearby absorbers. 



5.1. Uncertainty in the CfA2 Luminosity Function 

The uncertainty in the CfA2 luminosity function (MHG94) introduces an 
uncertainty in our selection function weighting (eq. |3-4[ ). Our results for (3 and 
/g a i are insensitive to the overall LF normalization, 0*; it cancels out of equations 
( [4-1] ) and ( [4-3D . We may therefore restrict our error analysis to uncertainty 
in the LF shape parameters, a and M*. To incorporate this LF uncertainty 
into our error bars on (3 and / ga i, we Monte Carlo resample the densities in 
equations ( |4-1| ) and ([4-3|) , varying the LF parameters according to their ranges 



from MHG94: M* = -18.8 ± 0.3, and a = -1.0 ± 0.2. The la bounds on the 
absorber distribution parameters increase only marginally from §|j: 

(3 : -0.02 ± 0.23 -> -0.02 ± 0.24; / gal : O.OOtJj 24 -> O.OOig' 25 . (5-1) 

MHG94 show that the respective LFs of CfA2 North and CfA2 South differ 
significantly (> 2a) with the LF of the Combined (CfA2 North+South) sample. 
This North/South LF difference is difficult to interpret; it may be caused by 
large-scale systematic errors in the Zwicky photometry or by real variations in 
the underlying LF across the survey volume. We investigate its effect on our 
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density distribution analysis by replacing the selection function derived from 
the Combined sample (0* = 0.04 (Mpc//i)~ 3 , M* = -18.8, a = -1.0) with a 
hybrid selection function using the LFs from CfA2 North (0* = 0.05 (Mpc//i) -3 , 
M* = -18.67, a = -1.03) and CfA2 South (0* = 0.02 (Mpc//i)~ 3 , M* = -18.93, 
a = —0.89) as appropriate. The hybrid selection function results in insignificant 
changes to the model parameters: 

(3 : -0.02 ± 0.23 -> -0.05 ± 0.26; / ga i : 0.00^' 24 -> O.OO^ 25 . (5-2) 

Based upon Equations ( |5-1] ) and ( |5-2| ) , we conclude that our density neighbor- 
hood technique is highly robust against uncertainty in the survey LF. 

5.2. Choice of Smoothing Length 

To determine the sensitivity of our results to the particular galaxy smoothing 
length, we repeat the analysis of §^j with smoothing kernels of 2.5h~ 1 Mpc and 
10h~ l Mpc. These values are near the functional limits on kernel size: broader 
kernels wash out the large-scale structure we are hoping to probe, and narrower 
kernels are subject to excessive shot noise. 

Figures |8] and |9|, respectively, show the (3 and / ga i likelihood curves for kernel 
smoothing lengths of 2.5h~ l Mpc (dashed), 5/i _1 Mpc (solid), and 10/i _1 Mpc 
(dotted). As might be expected, the low-contrast lOh^ 1 Mpc smoothing length 
results in poorer parameter constraints; the high contrast of the 2.5h~ 1 Mpc 
smoothing gives tighter constraints. These curves do not include shot noise cor- 
rections (§ |5.3j) , which have a larger broadening effect at the smaller smoothing 
length. We give the la confidence intervals on /3 and f ga \ for the three smooth- 
ing lengths in Table ^. For the range of smoothing length 2.5— 10/z — 1 Mpc, the 
maximum likelihood values for (3 and / ga i remain consistent with (randomly 
distributed absorbers) at the 1.2a level or less. Likewise, this range of smoothing 
lengths remains inconsistent with (3 = 1 (absorbers tracing structure) at the 3.8cr 
level or more. Thus our results are largely insensitive to the particular choice of 
smoothing kernel scale over a reasonable range. 

5.3. Shot Noise in the Density Estimator 

As a first estimate of the (n/n) shot noise, we construct a mock redshift survey 
drawn from the CfA2 LF. We place the mock survey galaxies at random locations 
in redshift space, with the underlying density (n/n) = 1 throughout. We then 
remove galaxies fainter than the CfA2 magnitude limit and use the remaining 
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sample to estimate the density contrast throughout the survey volume. By con- 
struction, the (n/n) shot noise throughout this homogeneous survey should only 
be a function of redshift. Figure ITDI shows the variance of the density estimate, 
binned in 2.5h~ l Mpc radial intervals out to 150/i _1 Mpc. Shot noise rapidly 
dominates the density estimate for cz ^ 12000 km s _1 , where the galaxies on 
average are separated by more than twice the 5h~ l Mpc smoothing length. Thus 
we restrict our analysis to Lya absorption systems with cz <^ 10500 km s _1 . The 
systematic decline in the median estimated density at random locations with 
cz ^ 10000 km s _1 has virtually no influence on our test population of random 
locations (§ |4.3.1 ) , drawn from the cz max <^ 10500 km s _1 radial intervals of Table 

Unlike our mock survey, CfA2 is not homogeneous. The true shot noise will 
be relatively larger in regions of low galaxy number density, and vice versa. We 
reproduce this structure-dependent shot noise by bootstrap resampling the entire 
redshift dataset. We recompute the likelihood (3 and / ga i curves (eqs. |4-1| and fi~3| ) 
for 1024 bootstrap realizations of the redshift dataset. The la bounds expand to 
-0.26 < (3 < 0.22, and 0.00 < / ga i < 0.25, as compared with the ranges in § |4~3| . 
This minor increase grows large if we attempt to include the next two nearest 
Lya systems, which are in the 3C 273 foreground at cz m 15000 km s _1 . The 
large uncertainty (±1 dex) in the density estimate at 15000 km s _1 overwhelms 
any benefit to be gained from the larger sample size. In §5"!4", we discuss the 
observational requirements to include the 3C 273 absorbers at z ^ 0.05 in our 
analysis. 



5.4. Improving the Parameter Constraints with a Larger Sample 

With 18 Lya absorbers, the distribution of surrounding galaxy density con- 
trasts places interesting constraints upon the nature of the low-redshift Lya for- 
est. A larger sample of absorber local densities should further shrink the error 
bars on the parameter constraints. Furthermore, local Lya absorber searches 
along additional LOS would improve the likelihood that the cumulative path- 
length probes a representative distribution of the galaxy number density within 
the survey. 

We quantify the improved constraints with a Monte Carlo simulation involv- 
ing 50 mock absorbers placed along random LOS with the same mean frequency 
as our current sample. We assign redshifts to the mock absorbers according to 
an a priori (3- or / ga i-model. The density contrast analysis recovers the seed val- 
ues (3 or / ga i with la error bars of only ~0.10. Were the nearby Lya absorbers 
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truly distributed randomly, for example, a distribution of density contrasts for 50 
absorbers should conclusively rule out a model with /3 — 1 at the 9a level. Assum- 
ing a local absorber frequency of ~ 1/2500 km s -1 detectable with HST/STIS, 
a sample of 50 absorber densities would require roughly a dozen AGN lines of 
sight within surveyed regions. 

The current limits to the sample size of local Lya absorber densities are: 
(1) the small number of AGN sightlines with identified low-redshift absorbers 
and (2) the limited depth of the complete redshift surveys for density estima- 
tion. As an example of the latter, there are ten more 3C 273 absorbers, with 
0.05 ^ z ^ 0.16, whose local densities cannot be estimated from a m Zw < 15.5 
redshift survey. Including all the additional 3C 273 absorbers would increase the 
sample size to 28 and almost double the cumulative pathlength currently probed 
by the various LOS in this study. Below we make a rough estimate of the ob- 
servational requirements for studying the large-scale density environment around 
these more distant absorbers in a manner consistent with our current sample. 
Our technique differs substantially from the approach suggested by Sarajedini, 
Green, & Jannuzi (1996) for the even more distant (0.2 < z < 0.4) absorbers 
detected with HST/FOS. 

To determine the (n/n) profile along the entire 3C 273 LOS without significant 
edge dampening or shot noise, we need to survey galaxies brighter than ~ 
out to ~ 2 smoothing lengths, or lOh^ 1 Mpc, from the LOS. At z — 0.05 this 
corresponds to B ~ 16.5 over ~50 square degrees, decreasing to B ~ 19.5 over 
~ 5.6 square degrees at z — 0.16. For comparison, Morris et al. (1993) have 
already surveyed ~3.5 square degrees around 3C 273 to B ~ 19, and Grogin et 
al. (1998) have recently completed a 2 h x 12° redshift survey around 3C 273 to 
the Zwicky catalog limit (m Zw = 15.7). 

6. Conclusions 

We develop a new diagnostic for exploring the relationship between Ly« 
absorption systems and the galaxy distribution: we compare the CDF of the 
smoothed galaxy number density at the absorber locations with the density CDF 
for test populations within the redshift survey. The underlying assumption is that 
if Lyct absorbers sample the large-scale structure revealed by redshift surveys, the 
distribution of smoothed galaxy number densities around absorber positions in 
redshift space should be consistent with densities around galaxies in the redshift 
survey. 

The procedure for evaluating the smoothed density field is straightforward, 
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requiring only the galaxy positions and redshifts, the survey selection function, 
and a choice of smoothing kernel. Moreover, this technique may be easily gen- 
eralized as long as the redshift survey dimensions are large compared to the 
smoothing kernel. For sensitivity to large-scale structure without excessive blur- 
ring of the survey density contrasts, we adopt a Gaussian smoothing kernel of 
width ^h^ 1 Mpc, which is approximately the galaxy-galaxy correlation length. 

The density CDF of 18 Lya absorbers within regions surveyed to mz w = 15.5 
is inconsistent with the comparable CDF for CfA2 galaxies. The K-S probability 
that the two CDFs are drawn from the same underlying distribution never exceeds 
2 x 10~ 4 for a variety of CfA2 sampling geometries and weighting schemes (§ §4.1| . 
[4.2j ). The low probability results from a much larger fraction of Lya absorbers 
than galaxies in underdense regions of the redshift survey. By contrast, a CDF 
of densities at randomly selected locations in the survey represents the same 
underlying density distribution as the Lya absorbers with a K-S probability of 
20%. 

We explore the likelihood that the nearby Lya absorbers are: (1) a single 
population distributed as the galaxy density raised to some power-law (3, or (2) 
a dual population of absorbers, with a fraction / ga j distributed like the galaxies 
((3 = 1) and the remainder distributed uniformly with respect to the galaxies 
((3 = 0). For the 18 absorbers in this study, the ±1<t likelihood interval of (3 is 
—0.02 ± 0.26. The error bars include the uncertainty in the density estimator 
from the survey sampling variance and the uncertain survey LF, both small 
contributions. The maximum likelihood value of f ga i is 0.0, with a la upper limit 
of 0.26. This value is consistent with the / ga i ~ 0.20 found by Mo & Morris 
(1994) for the 3C 273 absorbers alone, but at odds with the / ga i = 0.65 ± 0.18 
found by Lanzetta et al. (1995) for a sample of HST/FOS absorbers at higher 
column density. 

We conclude that the low redshift, low column density Lya clouds in this 
sample are not tracing the nearby large-scale structure marked by typical lumi- 
nous galaxies. The absorbers appear in a range of density environments similar 
to those around randomly chosen locations throughout the survey. Thus the 
Nhi 10 14 cm~ 2 clouds at z < 0.035 are apparently similar in their spatial 
distribution to the unclustered, low column density clouds seen at high z ( (5ar-| 
gent et al. 1980| ). The sharp decline in absorber metallicity at <^ 10 14 re- 
cently reported by Lu et al. (1998) suggests that the absorbers found nearby with 
HST/GHRS may be a distinct population from the HST/FOS absorber sample, 
which Lanzetta et al. (1995) argue is associated [f ga i ~ 0.65) with ~ 160/i -1 kpc 
galactic halos. 
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Although our sample is large enough to reject / ga i = 1, it does not exclude 
the conservative lower limit, / ga i = 0.35±0.10, of Lanzetta et al. (1995). In order 
to reduce the f3 and / ga i error bars to ±0.10, the precision necessary to exclude 
the Lanzetta et al. lower limit, Monte Carlo simulations indicate that a sample 
of ~ 50 nearby absorbers is required. The additional ~ 10 sightlines needed for a 
50-absorber sample would also increase the likelihood that the cumulative path- 
length intersects the full range of density environments in the nearby universe. 
With the progress of wide-field imaging and multi-object spectroscopy, redshift 
surveys around the HST Lya absorber sightlines using 4m-class telescopes should 
soon enable the density environment analysis presented here to encompass many 
more Lya absorbers out to z <^ 0.4. Investigation of absorber density CDFs for 
subsamples grouped by column density might then reveal a contrast in absorber 
spatial distribution akin to the metallicity falloff seen by Lu et al. (1998). 

We thank Susan Tokarz, Emilio Falco, and Michael Kurtz for assistance with 
the CfA2 survey database, and Jan Kleyna for useful discussions. This research 
is supported in part by the Smithsonian Institution. 
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Table 1: Targets of Local Lja Absorption Searches 



Target 


R.A. 


Dec. 


Redshift a 


Lya Sensitivity 3, 




(B1950) 


(B1950) 


(km s- 1 ) 


(km s- 1 ) 


Mrk 335 


00 h 03 m 45!2 


+ 19°55'29" 


7901 


1713-7901 


I Zw 1 


00 h 50 m 57!8 


+ 12°25'20" 


18490 


1660-10660 


Mrk 421 


11 h 01 m 40 s 6 


+38°28'43" 


9001 


1501-9001 


3C 273 


12 h 26 m 33!2 


+02°19'43" 


47347 


0-10500 b 


Mrk 501 


16 h 52 m lr!7 


+39°50'25" 


10301 


1709-10301 



a Galactocentric velocities. 

b Absorbers detected to z~ 0.12, but CfA2 too sparse for meaningful density estimation. 



Table 2: Local Lya Absorption Velocities and Equivalent Widths 



Background 


C^abs 


w x 


Background 


C^abs 


W x 


Source 


(km s _1 ) 


(mA) 


Source 


(km s _1 ) 


(mA) 


Mrk 335 


2183 


170 


3C 273 


2060 


240 




2503 


73 




6413 


27 




4483 


26 




7737 


32 




6493 


140 




8704 


120 


I Zw 1 


1777 


120 




9750 


74 




3021 


84 


Mrk 501 


4869 


154 




5290 


84 




6209 


36 


Mrk 421 


3047 


92 




7739 


48 


3C 273 


890 
1460 


371 
414 









Note. — Velocities are galactocentric. 3C 273 data from Morris ct al. (1991), Bahcall et al. (1991), and 
Weymann et al. (1995); remaining data from Shull, Stocke, & Penton (1996). 
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Table 3: Results for Various Smoothing Kernel Scales 

Param. Kernel Smoothing Length 

2.5 Mpc 5 Mpc 10 Mpc 

P 0.09 ±0.13 -0.02 ±0.23 -0.45 ± 0.38 

fa o.Q6±8;gg o.oo±g;gg o.oo±8:gg 



Note. - 



la (68.3%) confidence intervals, without corrections for shot noise or LF uncertainty. 
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Figure Captions 



Fig. 1. — Redshift survey slice and density profile along 3C 273 line of sight. 



See §3.2 for description. 



Fig. 2. — Redshift survey slice and density profile along Mrk 335 line of sight. 



See § |3T2| for description. 



Fig. 3. — Redshift survey slice and density profile along Mrk 421 line of sight. 



See §3.2 for description. 



Fig. 4. — Redshift survey slice and density profile along Mrk 501 line of sight. 



See § |3T2| for description. 



Fig. 5. — Redshift survey slice and density profile along I Zw 1 line of sight. See 
§|3.2| for description. 



Fig. 6. — Cumulative distribution function of galaxy density contrast around 18 
nearby Lya absorbers (solid), and around all CfA2 galaxies. The CfA2 sample's 
CDF is shown with (dotted curve) and without (dashed curve) selection-function 
weighting (cf. § |4.1|) . 



Fig. 7. — Cumulative distribution functions of galaxy density contrast for var- 
ious samples probing sightlines through the CfA survey, including: 18 nearby 
Lyct absorber locations from five AGN sightlines (solid); CfA2 galaxies within 
2.5h~ 1 Mpc of random sightlines, with (dotted) and without (short dash) 
selection-function weighting (cf. § |4.2|) ; and random redshifts from random sight- 
lines, with the same mean redshift spacing as the Lya sample (long dash). 

Fig. 8.— Likelihood curves for models with Lya absorbers distributed as a 
power law (3 of the galaxy density contrast (cf. § [4.3. 1|) . Results shown for galaxy 
smoothing lengths of 2.5h~ l Mpc (dotted), 5/i -1 Mpc (solid), and \0fi~ 1 Mpc 
(dashed). 

Fig. 9. — Likelihood curves for models with Lya absorbers distributed in two 
populations, with a fraction / ga i distributed as the CfA2 galaxies and the re- 
mainder distributed randomly (cf. § |4 . 3 . 2| ) . Results shown for galaxy smoothing 
lengths of 2.5/i~ 1 Mpc (dotted), hh' 1 Mpc (solid), and 10/z 1 Mpc (dashed). 

Fig. 10. — Variance in the 5/i -1 Mpc density contrast, {n/n), within a simulated 
CfA2 survey of uniform density throughout (dotted line). The median value 
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(crosses) and ±ler bounds are shown for the density contrast sampled at a large 
(N > 1000) number of random locations within successive 2.5/i _1 Mpc radial 
shells at distance D. Because of the growing systematic and statistical errors 
in the density estimator at cz ^ 12000 km s _1 , we limit our analysis to Lya 
absorbers with cz < 10500 km s _1 . 



FIGURE 1 




(u/u) 



FIGURE 2 




(u/u) 



FIGURE 3 




(u/u) 



FIGURE 4 




(u/u) 



FIGURE 5 



o 




(u/u) 



FIGURE 6 



c 

.2 0.8 

-*-> 

o 

c 

.2 0.6 

-l-> 

p 

. I — I 

• i — i 

Q 



0.4 



0) 

> 

. , — , 

-+-> 

cd 



5 0.2 







0.1 



"i i r 



/ 



-4- . 1--'"' I I I L 



J I I I I I I L 



(n/n) 



10 



FIGURE 7 




FIGURE 8 




-1 1 





FIGURE 9 




FIGURE 10 



i r 



z i- s- t § i i i -i -J f-fT c 



1.1 



I I I I I I 

50 100 



J I I L 



D [h- 1 Mpc] 



150 



